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Atomic Bose-Einstein condensates confined to a dual-ring trap support Josephson vortices as 
topologically stable defects in the relative phase. We propose a test of the scaling laws for defect 
formation by quenching a Bose gas to degeneracy in this geometry. Stochastic Gross-Pitaevskii 
simulations reveal a —1/4 power-law scaling of defect number with quench time for fast quenches, 
consistent with the Kibble-Zurek mechanism. Slow quenches show stronger quench-time dependence 
that is explained by the stability properties of Josephson vortices, revealing the boundary of the 
Kibble-Zurek regime. Interference of the two atomic fields enables clear long-time measurement of 
stable defects, and a direct test of the Kibble-Zurek mechanism in Bose-Einstein condensation. 



A possible mechanism for the formation of domain 
structures in the early universe was proposed by Kib- 
ble [l|. In this original paper, it was argued that the 
Universe cooled down after the hot Big Bang event and 
subsequently passed through a symmetry breaking phase 
transition with critical temperature Tc. Causally un- 
connected spatial domains settling into different vacua 
would lead to the formation of defects like domain walls, 
monopoles, strings, textures, etc [2]. Due to thermal fluc- 
tuations thwarting the emerging order, it was postulated 
that the number of defects eventually settled at the so- 
called Ginzburg temperature Tq < Tc. 

Later Zurek jsf put forward an alternative argument 
that focused on the nonequilibrium aspect of the phase 
transition. The density of the defects is determined at 
the critical temperature instead and its number is scaled 
with the quench rate. The scaling exponent depends 
on the critical exponents of the underlying phase transi- 
tion. This scenario, known as the Kibble-Zurek mecha- 
nism (KZM) should equally apply to condensed matter 
phase transitions where it is accessible to laboratory ex- 
periments [^]. The KZM proved to be robust and was 
verified by a number of recent experiments on supercon- 
ductors jsj-lsl and theoretical research on Bose-Einstein 
condensates (BEC) [l, [III, EH, EH- It also extends to 




quantum phase transitions 

In the spirit of Kibble's original argument one might 
expect the KZM to fail in the limit of slow quenches 
where the time scale of other processes occurring in the 
system dominates over the quench time. Deviations from 
KZM predictions were observed in ^He experiments [l6| 
but interpretation was controversial [l^ and a manifes- 
tation of the Ginzburg temperature was ruled out in 
Ref. [l8|. So far a transition between a regime of KZM 
scaling and its breakdown has not been studied system- 
atically. 
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FIG. 1. (Color online) Schematic of the two linearly coupled 
BECs. The isosurface shows the equilibrated condensate den- 
sity profile for a double-ring trap configuration and the color 
shows a phase profile with three Josephson vortices resulting 
from a quench at fixed tunneling energy J. The solid line 
on the left visualizes the trapping potential in the vertical di- 
rection, which forms a double well. The interference pattern 
of the two atomic fields on the bottom shows clear evidence 
of the Josephson vortices which are characterized by the low 
density regions. 



In this Letter we investigate the robustness of the KZM 
in a model system where departure from KZM scaling 
can be understood in detail because the defects are eas- 
ily quantified and are stable at the end of the quench. 
This avoids the difficulty of counting the decaying popu- 
lation of defects [ll|, [20| or their remnants [l9| . To this 
end, we study two linearly coupled quasi-one-dimensional 
(quasi- ID) atomic Bose gases in the ring configuration, 
as in Fig. [TJ A quench through the Bose-Einstein con- 
densation phase transition can generate Josephson vor- 
tices ( JVs) confined between the two BECs [2l|, |22| . We 
show that the number of these defects obeys the Kibble- 
Zurek scaling law for fast quenches. On the contrary, 
for slow quenches, the predicted behavior deviates sub- 
stantially and we observe a much stronger quench-time 
dependence than expected for critical phenomena in our 
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simulations. This is due to decay processes occurring 
before topological stability is established, in analogy to 
Kibble's arguments. 

We consider two linearly coupled Bose gas in the quasi- 
ID ring geometry. The trapping potential can be real- 
ized by crossing a vertical Gaussian-Laguerre laser beam 
(LGj) and two horizontal sheet beams [2^ to form an op- 
tical dipole trap or with rf-dressing on an atom chip [23| . 
Along the z-axis the trapping potential can be treated as 
a double- well potential as shown on the left of Fig. [TJ As- 
suming tight confinement, the transverse degrees of free- 
dom can be eliminated. The resulting coupled quasi- ID 
Gross-Pitaevskii equations for the BEG order parameter 
ipi and 7/^2 in each ring assume the dimensionless form 



(1) 



where Cj 



-\dxx and j = 1,2. Length (x), 

time (t), and energy (/i, J) are given in the harmonic os- 
cillator units of ah = ^^h/muo^ th = and Eh = ^cj, 
respectively, where m is the atomic mass and uj is the 
trapping frequency in the transverse direction. Dimen- 
sionless g = 2a /ah is the non-linear interaction strength 
with the 5- wave scattering length a; fi is the chemical 
potential, and J is the tunneling energy. 

Equation ([T]) supports topological and non-topological 
defects in the form of the JV and the dark soliton, re- 
spectively. 



^^1,2 = Vl + ^ tanh (px) ± iBsech (px) , 



(2) 



where u = J/fi^ and the scaling x = ^/Jlx and i/jj = 
y/jjjg'ipj has been applied. Both the dark soliton with 
B = and p = \J\^v and the JV with p = 2y/v 
and B = Vl - Zv for v < 1/3 are excitations above 
the vacuum where = ^p2 = const and are localized 
on the length scale ah{^/~jjip)~^ [Uj. The dark soliton, 
where both components have identical profiles, is non- 
topological because it is continuously connected to the 
vacuum by a family of moving "grey" solitons with de- 
creasing energy [25| . Although they may be present tran- 
siently during a quench through the BEG phase transi- 
tion, dark solitons will thus not survive the final stage 
of cooling. Furthermore, for v < 1/3, dark solitons are 
dynamically unstable with respect to the decay into JVs, 
which have lower energy [26|. The stability properties of 
the JV, on the other hand, depend on the dimensionless 
parameter u and may change during the quench. The 
JV bifurcates from the dark soliton at = 1/3 as a time- 
reversal symmetry broken state (vortex and anti- vortex) 
with a characteristic phase winding of 27r around a point 
located between the two rings (see Fig. [T]), and only ex- 
ists for smaller values of u. From numerical simulations 
it is known that JVs can move with respect to the back- 
ground BEG, although explicit solutions are unknown. 
For 1/5 < ly < 1/3 perturbative arguments indicate that 
the JV is energetically unstable [26|. For u < 1/b where 



the JV resembles the Sine-Gordon soliton [2l|, |22|, the 
stationary solution (|2]) is a metastable local energy min- 
imum, since the energy increases with velocity. Thus, 
at sufficiently small JVs are topologically stable, en- 
abling experimental tests of the KZM scaling by counting 
the number of JVs at the end of the quench in a dual- 
ring BEG. The defects would be immediately evident by 
interference of the two atomic fields via expansion imag- 
ing. The situation is strikingly different to a single one- 
dimensional BEG where the KZM scaling law was pre- 
dicted to govern a transient population of dark solitons 
that finally decays, which makes experimental detection 
more difficult MJJ. 

The nonequilibrium dynamics during the thermal 
quenches can be described by the coupled stochastic 
Gross-Pitaevskii equations [13, [lH : 

dt/jj = (i + F) [(/i(t) - Cj)tljj + Jil^s-j] dt + dWj, (3) 

where F is the growth rate, and fi is the chem- 
ical potential of the thermal cloud. dWj is the 
white noise associated with the growth process, which 
is characterized by the ffuctuation-dissipation relation 
{dW^ (x, t) dWk (x^ t)) = 2rTSjkS {x - x') dt, where T 
is the temperature in units of hw/kB- At the mean 
field equilibrium level the phase transition is described 
by the ground state of the energy % = j dx ^\dx'4^i\^ ^ 
\\dx'4^2\^ + y{'4^ii'4^2)^ where we seek the minimum of 
the potential V(V^i,V^2) = Ej=i,2 IV^i P [f l^^il^ " " 

+ V^2V^i] for J > 0. The symmetry V(V^i,V^2) 
y{i^2i'^i) imposes a common amplitude for the ground 
state fields. Taking = ^/ne^^^ , ^^2 = \fne^^'^ , and 
A = 6^1 — ^2, the minimum of F = gv? — 2\in — 2Jn cos A 
occurs at A = 0, n = {ji^ J) / g, ioi ji > — J. At the crit- 
ical point ji = — J the minimum is independent of A and 
each field can acquire a [/(I) symmetry-breaking phase. 

The dynamical transition to the broken symmetry 
phase is simulated via Eq. (j3j) with time-dependent chem- 
ical potential 



(4) 



where tq is the quench time. The quench starts from a 
thermal gas with chemical potential — /io < 0, and ceases 
in the Bose-condensed phase at /io > 0- Due to the 
coupling between the two rings jl{t) = t / tq -\- J acts as 
the effective chemical potential in mean field theory; the 
precise location of the transition in a dynamical quench 
must be determined numerically. To count the number 
of JVs during the quench, we evaluate the total phase 
defined as 6^^ = ^ dOi for the i-th ring. Gonsequently 
the total number of JVs can be quantified by calculating 

^"7^ = / 1^(^1-^2)1/2^. 

Kibble-Zurek theory applied to the BEG phase transi- 
tion gives the relaxation time and healing length close to 
the critical point as 



-1/2 



(5) 
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FIG. 2. (Color online) Scaling of the total number of JVs 
(Njv) with respect to tq at J = 5 in (a) and J = 25 
in (b) averaged from 500 trajectories of Eq. (|3]). The er- 
ror bars indicate the standard deviations. The red lines 
show the best power-law fit for fast quenches with exponents 
-0.2523 ± 0.0128 in (a) and -0.2456 ± 0.0131 in (b). These 
fit very well with the KZM prediction of -1/4. The dashed 
lines indicate the critical quench time Tq^^ of Eq. (p^ for the 
breakdown of the KZM scaling law. 

where and tq depend on the microscopic details of the 
system. Using Eq. dH and following the Kibble- Zurek 
scenario [29], we obtain the typical size of the domains 
after the quench 

where for our system tq = F"^. When /i (t) exceeds 
— J, localized phase-domains start to grow in each ring. 
Typically a piece of the (anti-) vortex will fall within a 
^-sized domain and the phase in each domain is chosen 
independently and randomly. Therefore, for small J, in 
a ring with circumference C, the total number of JVs is 
estimated to be 

(iVjy)~C/e = Ceo-^(^) (7) 

and thus obeys the —1/4 power-law scaling with quench 
time. The number of JVs thus shows a stronger quench- 
time dependence than the winding number of a single- 

2/8 

ring BEC, which was predicted to scale with Tq ' (15l] . 

We consider a gas of ^^Rb atoms with a transverse con- 
fining frequency of a; = 27r x 200Hz and corresponding 



nonlinear interaction strength ^ = 0.05. To study the 
quench dynamics we numerically integrate Eq. (j3j) with 
circumference C = 30 and T = 10~^. All these param- 
eters are realistic, e.g. with the set-up of Ref. [30|. The 
scaling in Eq. ([7]) is verified by averaging N jy over 500 
trajectories for two tunneling energies, J = 5 and 25. 
The quenches vary the chemical potential from — /io to 
+/io, with /io sufficiently large that the resulting defect 
number is independent of it. As shown in Fig. [21 the sim- 
ulation results of fast quenches compare favorably with 
the KZM prediction, however the number of JVs deviates 
from the KZM scaling for slow quenches. 

The stability of a JV depends on conditions that 
change during the quench. According to the KZM, two 
different regimes exist: For early as well as late times dur- 
ing the quench, relaxation is efficient and fluctuations in 
the Bose gas follow the changing chemical potential adi- 
abatically. However, when the diverging relaxation time 
r of Eq. ([5]) exceeds the time scale of the quench /i//i, 
fluctuations transiently freeze out and the system enters 
the impulse regime. This occurs when 

T{fi{i)) = = i, (8) 

giving the freeze-out timescale t = ^Jtqtq. At the follow- 
ing impulse-adiabatic transition the frozen fluctuations 
are imprinted onto the forming BEC. We thus expect that 
the stability properties of formed defects at this transi- 
tion point determine their survival during the adiabatic 
phase of the quench. 

A convenient observable to detect the impulse- 
adiabatic transition is the particle number, which is 
shown in Fig. [3l While the particle number is small 
in the impulse regime, it follows the dashed linear time- 
dependence in the adiabatic regime. Figure [3] shows clear 
evidence of the impulse-adiabatic transition happening at 
ji = fi/rq^ with a value of / ^ 4.7 read off the numer- 
ical data; the value of / appears to depend only weakly 
on the details of the system, including the parameters J 
and F, consistent with the theoretical argument that ji is 
relevant for the quench dynamics. Similar behavior was 
observed for a single ring BEC Therefore we can 

predict the chemical potential at the impulse-adiabatic 
transition as: 

fl = f^{ft-TQj)=f^-J. (9) 

We denote the critical ratio of tunneling to chemical po- 
tential for JV stability by Vc = J/fJ^c^ where fic is the 
stabilizing chemical potential for given J. As shown in 
Fig. IHa), the defects are frozen in until fi > jic for fast 
quenches ensuring the topological protection of JVs and 
hence the KZM signature. However, for slow quenches 
the impulse regime terminates earlier with jl < fic-, which 
causes the decay of the JVs in the shaded region in 
Fig. m until the topological stability of JV is established 
at ii{t) = lie Although the the stability properties of 
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FIG. 3. (Color online) Particle number of component '0i as a 
function of time for J = 5 in (a) and J = 25 in (b). The in- 
set enlargements have the vertical scale magnified by a factor 
of 10 and reveal details for slow quenches. The color-coded 
labels show tq. Quenches with vastly different tq show a 
knee structure characteristic of the impulse-adiabatic transi- 
tion with a rapid particle number increase around jiTQ/t = 4.7 
for both J = 5 and 25. 



the static JV at equilibrium are known, the critical Uc 
for a moving JV at finite temperatures remains unclear. 
Nevertheless we can estimate from the numerical sim- 
ulations, as the point where the JV becomes unstable. 
From Eq. ([9]) we obtain the criterion for obtaining stable 
JVs 



< Tq 



(10) 



The value Vc — 0.0813 is obtained from the data for J = 
5. With this value at hand, we can predict Tq'^^ = 0.02 
for J = 25, as shown by the vertical dashed line plotted 
in Fig. [2jb). This prediction agrees with the numeri- 
cal data very well. The critical quench time depends on 
the growth rate through tq and we have also verified the 
prediction of Eq. (p!Q|) at different growth rates. In the 
slow quench regime the defect number varies much more 
strongly with quench time than expected from the KZM 
as is seen by the negative gradient of large magnitude 
seen in Fig.O Since the gradient is clearly much stronger 
than the —1/4 expected from the standard argument be- 
hind Eq. ([7]), we do not expect to enter a new regime of 
power-law scaling for slow quenches. It should be noted 
that slow quenches show the same knee structure char- 
acteristic of the impulse-adiabatic transition as the fast 
quenches that lead to KZ scaling (Fig. [3]). We have also 
verified that the defect numbers for slow quenches con- 
tinue to satisfy the KZM scaling (solid line in Fig. 2) 



FIG. 4. (Color online) Schematic plots of the relaxation time 
r = ^o/lAl vs. chemical potential. During the quench, the 
system moves from left to right along the horizontal axis, jl 
is the chemical potential at the impulse-adiabatic transition 
and lie is the critical potential for stabilizing JVs, which are 
stable to the right of the dashed line. Panel (a) shows the fast 
quench where the KZM scaling is not affected by the stability 
of JV. Panel (b) represents a slow quench where the resulting 
number of defects is affected by decay processes taking place 
in the shaded area. 



when counted immediately after the impulse-adiabatic 
transition at fi. This supports our argument that the re- 
duced defect number is due to thermal decay processes 
happening after the transition, and that spontaneous de- 
fect formation is not affected by thermal fluctuations dur- 
ing freeze-out. 

For JVs, the regime for slow quenches is similar to 
the original idea put forward by Kibble [l|, where ther- 
mal fluctuations may still have enough time to destroy 
the emerging order before the decreasing temperature 
approaches the Ginzburg temperature. We observe a 
Ginzburg-like regime where thermal effects destroy the 
pattern of symmetry breaking inherited from criticality 
during the interval (i < fi < jic^ shown in the shaded 
region of Fig. IH^b). 

The absence of any striking, clear cut evidence of cos- 
mological nature and the difficulty in observing the KZM 
in condensed-matter experiments despite its dramatic 
predictions is usually not attributed to the failure of the 
mechanism but may be explained by the eventual decay 
of defects in the post-quench era [lil,|li,[20|. This diffi- 
culty is circumvented if the formed defects are manifestly 
protected from post-quench decay, e.g. by topology. The 
defects observed in the successful experiments of Refs. 
8] have this property as do the JVs that are the subject 
of this Letter. While cosmological defects protected by 
topology may still survive in dark matter or dark energy, 
their detection is difficult (3l| . 

Our work paves the way for a direct test of KZM in 
the Bose-Einstein condensation phase transition, by elim- 
inating post-quench decay of defects. The quench of ji 
could be supplanted by a controlled sweep of J, providing 
another knob for varying the quench time. 
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